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The recently derived long-range two-pion exchange (TPE) contributions to the nuclear current 
operator which appear at next-to-leading order (NLO) of the chiral expansion are used to describe 
electromagnetic processes. We study their role in the photodisintegration of and ^He and compare 
^3 ■ our predictions with experimental data. The bound and scattering states are calculated using five 

' different parametrizations of the chiral next-to-next-to-leading order (N^LO) nucleon-nucleon (NN) 

, potential which allows us to estimate the theoretical uncertainty at a given order in the chiral 

expansion. For some observables the results are very close to the predictions based on the AVI 8 
NN potential and the current operator (partly) consistent with this force. In the most cases, the 
addition of long-range TPE currents improved the description of the experimental data. 
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INTRODUCTION 



Chiral Effective Field Theory (ChEFT) provides a systematic and model-independent framework to analyze hadron 
structure and dynamics in harmony with the spontaneously broken approximate chiral symmetry of QCD. This 
approach is a powerful tool for the derivation of the nuclear forces. Exchange vector and axial currents in nuclei 
have also been studied in the framework of ChEFT. Since the pioneering work of Park et al. [l|, heavy-baryon 
chiral perturbation theory has been applied to derive exchange axial and vector currents for small values of the 
photon momentum. These calculations have been carried out in time-ordered perturbation theory. The resulting 
■ exchange vector currents have been, in particular, applied to analyze radiative neutron-proton capture within a 
, hybrid approach [2|. 

ChEFT has also been used to study the electromagnetic properties of the deuteron [3, 4] and ^He [H, Q. One of 
I the fundamental processes observed for the deuteron is the photodisintegration reaction. It has been a subject of 
intense experimental and theoretical research for several decades (se e Refs |3, HI). Also photodisintegration of ^He has 
been studied experimentally and theoretically for a long time [o'-lT] . Photodisintegration observables provide a good 
tool for studying the contributions from meson exchange currents (MEC) to the nuclear current operator. This is 
because the charge density operator, which often dominates low-energy electrodisintegration and is mostly given by the 
. single nucleon current, does not play any role in this reaction. An ongoing interest in low-energy photodisintegration 
J-H ■ reactions, especially in view of planned experiments, provides a strong motivation to apply the framework of chiral 
- - -' effective field theory. This approach relies on the approximate spontaneously broken chiral symmetry of QCD. It 
allows for a systematic derivation of the nuclear Hamiltonian and the corresponding electromagnetic current operator 
from the underlying effective Lagrangian for pions and nucleons via the chiral expansion, i.e. a simultaneous expansion 
in soft momenta of external particles and about the chiral limit. For more details on the application of ChEFT to 
nuclear forces and currents the reader is referred to recent review articles [12^, ^13] and references therein. 

In the two- and three- nucleon systems, the leading contributions to the exchange current originate from one-pion 
exchange which are well known. The 2N current operator at the leading loop order in the chiral expansion has been 
worked out by Pastore et al. [I^,[l5j based on time-ordered perturbation theory. Independently, the two-pion exchange 
2N current operator has been derived in Ref. [l6| using the method of unitary transformation. The resulting current 
operator is consistent with the corresponding chiral two- nucleon potential [12] obtained within the same scheme. 
In the present work, for the first time we explore the effects of the leading two-pion exchange 2N operator [l6| 
in the photodisintegration reactions of and ^He. We, however, emphasize that the presented calculations are 
not yet complete. In particular, the corresponding expressions for the one-pion exchange at NLO and short-range 
contributions to the current operator within the method of unitary transformation are not yet available. Our main 
goal in the present work is to explore the sensitivity of various observables in the deuteron and ^He photodisintegration 
to the two-pion exchange current rather than to provide a complete description of these reactions within the ChEFT 
framework. 
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Our manuscript is organized as follows. In section II the formalism which we use to describe selected 2N electro- 
magnetic reactions is presented. The results for the photodisintegration of the are discussed and compared with 
the experimental data in section III. The extension to the 3N system is briefly described in section IV and the results 
obtained for photodisintegration of ^He are presented in Section V. Finally, section VI contains the summary and 
conclusions. 



II. FORMALISM 



The general form of the nuclear matrix element for electromagnetic disintegration reactions in the 2N system is 
represented by 

m^{^Z,,\J>^m^lZnd), (2.1) 

where the proton- neutron scattering state l^scatt) deuteron bound state Inbound) obtained using NN 

potential. The current operator J^{Q) acts between the internal initial and final 2N states. We employ the solution 
of the Lippmann-Schwinger equation, t = V2N + tGoV2N^ in order to express as 

m = {po I (1 + tGo) J^{Q) I *^^u„d), (2.2) 

where Gq is the free 2N propagator, t is the NN t-matrix and |po) is the eigenstate of the relative proton- neutron 
momentum. Since all observables can be computed from ^ the description of the electromagnetic reactions requires 
the knowledge of the consistent potential and electromagnetic current. The NN potential based on ChEFT is currently 
available up to next-to-next-to-next-to-leading order in the chiral expansion [12, 13]. As already pointed out, in this 
paper we focus on the long-range TPE contributions to the current operator, which appear at NLO. However, in order 
to avoid the theoretical error from using the less accurate NLO NN potential, all calculations are made using the 
N^LO potential. At this order, the NN potential V2N is built from the one-pion exchange (OPE), Fi^^, and two-pion 
exchange (TPE), ¥27^^ contributions as well as various contact interactions (cont) [12| 

V2N = Vi^ + V2^ + Kont . (2.3) 

The effective current operator for the 2N system is a sum of the single- nucleon operators J^{i)^i = 1,2 and 
two-nucleon operators of different type (J^(l,2)) 

= J^(l) + J^(2) + J^(l, 2) , (2.4) 

where 



J^(l,2) = J(;(l,2) + J2^J1,2) + J,^„„t(l,2). (2.5) 

The expressions for the single-nucleon and the leading OPE currents Jf^(l,2) are well established, see e .g. |17| . The 
results for the leading two-pion exchange contributions used in the present work are available in Refs. [IJ^IlSl- We 
emphasize that the resulting two-pion exchange current is parameter-free. The expressions for the OPE and contact 
currents at the leading loop level have been recently worked out within time-ordered perturbation theory [|l5']. Work 
on the derivation of these contributions using the method of unitary transformation is still in progress. 

The 2N four-current operator J^(l,2) = ( J^(l, 2), J(l, 2)) can be decomposed according to its isospin and spin- 
momentum structure and quite generally written in the form [l6|, [l8| 

J^l,2) = ^^/f (gl,g2)T,0f , (2.6) 

11^1(3=1 
5 24 

^"^(^'2) = ^^/,^(gi,?2)T,d;3, (2.7) 

11^1/3=1 

where qi = p[ — pi is the momentum transferred to nucleon i, is the 2N isospin operator, 0| and O/3 are the 
(momentum dependent) spin operators in the 2N space, f^^ and are scalar functions. The explicit form of the 

scalar functions and the operator basis for and O/3 can be found in Ref. [16]. 

In this paper, we concentrate on a treatment of the long-range TPE contributions to the 2N current operator derived 
m Ref. 16 . The expressions for the functions f^^{qi, 02) and /^(^i, 02) entering the TPE current and charge density 
operators in Eqs. (|2.6p . (|2.7p are rather complicated and contain the standard loop functions and the three-point 
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functions in a form suitable for numerical calculations [l6|. Due to their isospin structure, not all combinations of 
(|2.6p and (|2.7p contribute to photodisintegration of the deuteron. The non- vanishing contributions emerge from 

10 

J2,(l, 2) = ^ f^{qiA2) (ri - T2)s Op + fi{quq2) (n x 7^2)3 62 , (2.8) 

(3=3 

where (...)3 denotes the third cartesian component of the vector. We work in momentum space and apply the standard 
partial wave decomposition of the 2N potential, see e.g. Ref. [19] for more details. Our calculations are performed 
using a complete set of 2N states 

\pa) = \p{ls)jmj)\tmt) (2.9) 

where p is the magnitude of the relative momentum, 5, j and mj are the orbital angular momentum, spin, total 
angular momentum and its projection on the quantisation axis i, respectively. The isospin quantum numbers of the 
two-nucleon system are denoted by t and m^. 

The TPE current operator needs to be expressed in the same partial wave basis. To this end, we first prepare all 
spin and isospin matrix elements using Mathematica© and then calculate the resulting four-fold angular integrals 

(pV|J2,(l,2)|pa) = {p\l's')fmjr^t'mt'\Jr,p\p{ls)jmj',tmt) (2.10) 
^ 0(r,'/;™„..»..-»..,m.)57,,,(p')C(i,i;m„m,-™„m,)y„„,(« 

mi 

X {01^02) {t'mt'\Tr^\tmt) {s' rrij' - mi'\di3\smj - mi) , 

numerically. Here, C{lsj;mi^mj — mi^rrij) denote the Clebsch-Gordon coefficient and YimiiP) are the spherical 
harmonics. Such an approach has been described in Ref. [20|. In order to calculate the four- fold integrals in Eq. (|2.10p 
for the whole grids of p and p^ points and all non- vanishing (a, a\ ruj) combinations we used the parallel supercomputer 
IBM Blue Gene/P of the Jiilich Supercomputing Centre (JSC). 



III. RESULTS FOR PHOTODISINTEGRATION OF THE DEUTERON 



We now discuss the results for the deuteron photodisintegration process for the unpolarized cross section and 
selected polarization observables. The results for the differential cross section, the photon analyzing power and 
outgoing proton polarization at the photon laboratory energies of = 10, 30 and 60 MeV are shown in Fig. [H The 
bands reflect the uncertainty due to the variation of the two cut-off parameters A and A that appear in the chiral 
potential. While the ffist cut-off parameter A appears in the regulator function for the Lippmann-Schwinger equation, 
the second parameter A enters the spectral function regularization (SFR) and denotes the ultraviolet cut-off value 
in the mass spectrum of the two-pion exchange potential. Following Ref. [l2[, the cut-off values are varied between 
450 and 600 MeV for A and between 500 and 700 MeV for A. It is important to emphasize that the resulting bands 
probably overestimate the theoretical uncertainty that can be expected in a complete calculation at this order in the 
chiral expansion. This is because we have not yet included the corresponding short-range contributions to the nuclear 
current, which are expected to absorb the large part of the cut-off dependence. Thus, the interpretation of the band 
width in the obtained results in terms of the theoretical uncertainty should be taken with care. Different bands shown 
in Fig. [1] describe the contributions from the different parts of the 2N current: single- nucleon current (light band), 
one-pion exchange contribution (hatched band) and the long-range TPE contributions (dark band). As a reference, we 
also show the results based on the phenomenological AV18 potential [21] and the corresponding current model j22l,l23|. 
Notice that the results obtained solely from the single nucleon currents and by adding the OPE contributions do not 
describe the data well and differ significantly from the reference AV18 predictions. An explicit inclusion of the TPE 
contributions yields an improved description of the experimental data which turns out to be in agreement with the 
AVI 8 predictions. The OPE predictions give the cross section and photon analyzing power values lower than the 
AV18 results. The bands including TPE currents are broad, however, for E.y = 30 and 60 MeV they give reasonable 
description of the experimental data. In the case of photon analyzing powers, the best agreement between all models 
is obtained at lowest energy. For energies — 30 and 60 MeV, calculations including the TPE currents yield even 
better agreement with experimental data than the AV18 results. It remains to be seen whether this conclusion will 
still hold after including the short-range and the subleading OPE currents. It is further important to emphasize that 
the single nucleon current alone is insufficient to describe the data for this observable. Thus, it is necessary to include 
higher order electromagnetic currents. In the case of the outgoing proton polarization, we observe a smaller sensitivity 
to the TPE currents and a good agreement between traditional framework (AVI 8) and chiral results at all energies 
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FIG. 1: (color online) The results for the unpolarized cross section, the photon analyzing power and outgoing proton polarization 
in the deuteron photodisintegration process at the photon laboratory energies of 10, 30, 60 MeV, displayed as functions 
of the proton emission angle. The solid black line refers to the standard calculation based on the AV18 potential, the light 
(blue) band covers results obtained with the single-nucleon current only, the hatched band represents the predictions based on 
the single-nucleon and OPE parts and the dark (pink) band includes, in addition, the contributions of the TPE current. The 
experimental data are from Ying et al. [5^. 



considered. The larger sensitivity to the details of the exchange currents is only observed at forward and backward 
outgoing proton angles. 

We have also calculated the deuteron tensor analyzing powers as a function of the proton emission angle for two - 
energy bins. Here, we focus on a comparison of our calculations with recent experimental tensor analyzing powers T2q 
for low energies from [24] and do not show the results for vector analyzing power iTu as there exist no experimental 
data for this observable. In order to be able to compare the theoretical calculations with the data from [2^, our 
predictions for the exclusive observables have been integrated over the relevant intervals of the initial photon energy 
and angular regions. In Fig. [2] the results for the angular distribution at the bin energies of = 2b — 45 MeV and 
E^ = 45 — 70 MeV together with experimental data are presented. For all deuteron tensor analyzing powers one 
observes a rather good agreement between the AVI 8 potential prediction, chiral results and experimental data. The 
effects of the TPE contributions turn out to be very small. Also, no broadening of the bands with increasing photon 
energy is observed. All this suggests that the deuteron tensor analyzing powers are driven by the long-range parts of 
the current and are not sensitive to the short-range contributions. We further emphasize some disagreement with the 
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FIG. 2: (color online) Deuteron tensor analyzing powers vs. proton emission angle for two E^- energy bins. The upper row 
shows results for bin energy — 25-45 MeV. The lower row shows results for bin energy E^ — 45-70 MeV. The bands and 
lines have the same meaning as in Fig. [T] The experimental data are from Rachek et al. ^] . 




FIG. 3: (color online) Total cross section for photodisintegration of the deuteron as a function of photon energy beam. In the 
left panel results for the single nucleon current and OPE contribution are shown. In the right panel, the results obtained with 
an additional TPE currents are given. The experimental data are the same as in Ref. 0]. 



data for T21 at forward angles. In the future, it would be interesting to see whether these conclusions are affected by 
the inclusion of the subleading OPE terms in the current operator. 

Finally, we have also calculated the total cross section. This observable was extensively studied in many publications 
using various theoretical approaches, see e.g. |7J. In Fig. [3l the total cross section for ^ 2-80 MeV is presented. 
The experimental data are taken from Refs. [3, [8|. In the left panel, we see the results for the single nucleon and 
OPE current. For this particular case the width of the prediction band is negligible. The theoretical predictions agree 
rather well with the experimental data. The right panel in this figure shows that the effects of the TPE currents are 
clearly visible, especialy at higher photon energies. We also notice that the band width increases significantly once 
the TPE contributions are included. 

To conclude, we observe that for all considered observables the 2N current operator plays an important role and 
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a restriction to the single nucleon current operator leads to a strong disagreement with the data. The inclusion of 
the leading OPE current is absolutely necessary to achieve a decent description of the experimental data. The effects 
of the TPE current are clearly visible in the differential cross section and some polarization observables such as 
and Py. One also observes a rather good agreement between the results based on chiral EFT and the AVI 8 potential 
combined with the corresponding current operators. 



IV. TWO-PION EXCHANGE CURRENTS IN THE 3N SYSTEM 

For 3N reactions, we use the framework and its numerical implementation described in detail in In this work, 
we will only briefly introduce the key points focusing mainly on the current operator in the 3N system. The starting 
point is exactly the same as for the 2N reaction. We consider the general matrix element of the current operator 
between the 3N bound state. Inbound)' scattering state l^scatt) system 

A^''^(*siIttk^(Q)l*bound)- (4.1) 

The 3N bound state Inbound) obtained in the standard way from the appropriate Faddeev equation [26*]. The 
current operator J^{Q) acts effectively between the internal initial and final 3N states. These internal states are 
conventionally expressed in the momentum space in terms of two Jacobi momenta, p and q [19]. The momentum p 
describes a relative motion within a 2N subsystem (here we choose the subsystem consisting of nucleons 2 and 3). 
The momentum q describes the motion of the spectator nucleon (here nucleon 1) with respect to that 2N subsystem 

P= ^{p2-P3^ , q = ^{pi - ^{P2 . (4.2) 

We consider two types of 3N scattering states. In the first case two nucleons bound in the deuteron emerge with 
the accompanying third nucleon and the asymptotic motion of this unbound nucleon is described by the Jacobi 
momentum ^o- In the second case, we have three free nucleons in the final state and their asymptotic relative motions 
are represented by p and q. 

In order to calculate the crucial matrix elements given in ()4.ip . it is not necessary to solve the corresponding 
Faddeev equations directly for the 3N scattering states [11]. Instead, we solve a Faddeev- type equation for an auxiliary 
state \ U). In our calculations, we do not include the effects of 3N forces. Thus, the equation for | U) takes a simpler 
form 

I U) = tGo (1 + Pyj'^iQ) I Hound) + tGoP I U) . (4.3) 
The corresponding nuclear matrix elements are then given by 

N^, = I (1 + P) J^(Q) I f g^,„d) + I P I [/) , 

Kpn = (*o I (1 + P)J''iQ) I *g^„„d) + (*o I (1 + P) I [/) , (4.4) 

with NN t-matrix in the 3N space. Go the free 3N propagator and P the permutation operator P = P12P23 + ^13^23- 
Further, is the antisymmetrized product state containing the deuteron and the momentum state for the relative 
motion of the third nucleon. Finally, |<I>o) is the antisymmetrized product state describing the two relative motions 
among the three outgoing nucleons. For details about the solution of Eqs (|4.3p - (|4.4p . see Ref. [11]. It is important to 
mention that these equations are solved in the partial wave basis. In the following, we, therefore, briefly discuss the 
partial wave decomposition of the current operator which can generally be written in the form 

= J^(l) + J^(2) + J^(3) + J^(2, 3) + J^(3, 1) + J^(l, 2). (4.5) 

There are three pairs in the 3N system, but it is sufficient to include a contribution just from one pair exploiting 
the fully antisymmetric nature of the 3N states. The two- nucleon current operator J^(2,3) is defined according to 
Eq. (2.5). Compared to Eq. ()2.8p . in the case of the 3N system, we have additional contributions from the Ti isospin 
structure. Thus, the non- vanishing contributions emerge from 

10 10 

J242,3) = ^ f?{q2^q3)Tid^ + ^ f^{q2^q3)T2d^ + fi{q2^ q3)Tsd2 , (4.6) 

(3=3 (3=3 

where q2 and qs are the momentum transfers of nucleons 2 and 3, respectively. We utilize the so-called jl coupling 
scheme for the 3N basis states 

\pqa) = \pq {ls)j(^X^)l{jl)JM)\(t^)TMT) =\pqaj)\aT). (4.7) 
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Here, /, 5, j and t refer to the orbital angular momentum, spin, total angular momentum and isospin of the (2-3) 
subsystem, respectively. The angular momentum of nucleon 1 is coupled with its spin ^ to the total angular momentum 
/. Finally, the subsystem total angular momentum j is coupled with / to give the total 3N angular momentum J 
with the projection M. Similar coupling in the isospin space leads to the total 3N isospin T with the corresponding 
magnetic quantum number Mr- 

Analogously to the procedure described in section II, we compute the general matrix element of the 2N current in 
the 3N basis 



'^27r(2,3) pqa^ = {p' q' aj>\d^ f^{q2,q^)\pqaj) {aT'\Trj\aT) 

= ^ C{j'l'J']mj>,M'-mj>,M')C{jIJ]mj,M-mj,M) 

rrij , rrijf 

X hsip'^P^Q; (l's')fmj',{ls)jmj) 

X h [q'. Q; I'M' - my, (^A0 IM - {ar' |T^| ar) (4.8) 



with 



h3{p\P:Q;{l's')j'mj^,{ls)jmj) = J dp' J dp^^C{l' s' j' ]mi' ,mj' - ,mj')Yi1 



(4.9) 



rriif mi 



C{lsj;mi,mj - mi,mj)Yi^rni{p) f^{q2,qz) {s' mj> -mv\dp\smj - mi) 



and 



h ?', ?, «i A'- I'M' -mfA\-\IM- m, 



jmtf-'-i,') 



where we have introduced 



(9) 



-v — m 
2 



(4.10) 



(4.11) 



For /23, we recognize the same type of a matrix element we dealt with in the 2N space. Now, however, the isospin part 
is separated out. This is because there are much more isospin combinations in the 3N system compared to processes 
on the deuteron which has the total isospin zero. This separation allows us, in particular, to calculate once and 
use it both for the reaction on ^He and ^H. For the numerical implementation it is, however, still important to use 
the properties of the matrix elements (q^t' I^t^I <^t) (r^ = 1,2,3) in order to reduce the number of necessary four- fold 
integrals in /23, even if the isospin dependence is now treated separately. Below we give the matrix elements of the 
three isospin operators in the 3N isospin space. As already mentioned, we assume that the operators act on the 3N 
bound state (^He or ^H), which has the total isospin T = 1/2. It is then straightforward to obtain 



t'\ 1 T'mT' 



1 \ 1 

*2 2™^ 



(t'i) T'mT'|(r(2)+r(3))3|(ii) \mT 



C(l,-,T';0,mr,mT 



I2^j{2t' + l){2t + l) 



1 t t' 

1 T-i/ 1 

2 ^ 2 



\t' t] 

(4.12) 



M 1 

*2 



i'i)TW' (r(2)-f(3))3 



C ( 1, -,T';0,mT,mT' 



1 \ 1 

*2 h"^r 



12 vW+l)(2i+T) 
(_l)t+t'+Hr' (^1 _ (_i)t+t' ^ ^ 



1 t t' 

1 T-T/ 1 

2 2 



(4.13) 
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FIG. 4: (color online) Differential cross section in the laboratory frame for ^He two-body photodisintegration at the photon 
laboratory energies = 12 MeV (left), = 20.5 MeV (middle) and E^ = 50 MeV (right). The band covers N^LO chiral 
predictions for different cut-off parameter values. The light (blue) band covers results obtained with the single- nucleon current. 
In the case of the hatched band, the current operator is taken as a sum of the single nucleons current and one-pion exchange 
current. The dark (pink) band covers N^LO chiral predictions for different cut-off parameter values and the current operator 
is taken as a sum of the single nucleons current, one-pion exchange current and two-pion exchange current. The solid line 
represents predictions obtained with the AVIS nucleon- nucleon potential and the related exchange currents [llj . 
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V. RESULTS FOR PHOTODISINTEGRATION OF ^HE 

We are now in the position to discuss our results for two- and three-body photodisintegration of ^He at three 
example photon laboratory energies E.y = 12, 20.5 and 50 MeV. The Coulomb force between two-protons in three- 
nucleon scattering states is not taken into account. The three- nucleon matrix elements N are obtained using the partial 
wave decomposition, with the total angular momentum of the three-nucleon system J < 15/2 and including all partial 
waves with the subsystem total angular momentum j < 3. The nuclear matrix elements TV are computed as using the 
formalism described in section IV. Given N ^ one can calculate cross sections and polarization observables which are 
expressed in terms of the nuclear matrix elements with different spin projections carried by the initial photon, the 
^He nucleus and the outgoing nucleons and/or deuteron. For more details we refer the reader to Refs. [H, |Tl|. In the 
following we just present our sample results for the chiral EFT approach. 

We begin with the exclusive unpolarized cross section for two-body breakup of ^He, d?(j/dVtdi where the final 
deuteron would be observed. It is depicted in Fig. |4] as a function of the deuteron scattering angle Od defined with 
respect to the initial photon direction at the photon laboratory energies E^ = 12, 20.5 and 50 MeV. We observe 
a similar behaviour as compared to the differential cross section in photodisintegration of the deuteron. The single 
nucleon current contribution yields significantly lower values as compared to the ones which include MFCs. The TPE 
bands overlap with the OPE bands and appear to be broader than the OPE bands. As expected, the bands become 
wider with increasing photon energy. 

Next, the results for a few polarization observables are shown in Fig. [5l We consider the photon {A2{0d)) and the 
^He {Ay^%Od)) analyzing powers as well as the spin correlation coefficients C2'^y^^{0d) and C^'^x^^i^d)- In the case 
of the photon analyzing power ^3^(7), the prediction bands for the single nucleon current give higher values than 
the other, more complete calculations, but the shape of the bands are always similar. The TPE bands are broader 
than OPE bands and overlap with them. For the ^He analyzing power ^^^(^He) and the spin correlation coefficients 
CxY and Cyx^ we observe that the results based on chiral EFT generate very broad prediction bands, especially at 
the highest energy considered. Interestingly, the results based on the single-nucleon current for these observables are 
completely different from the ones involving the MEG. This suggests that these observables are very sensitive to the 
details of the meson exchange currents, and their proper description will require the inclusion of the subleading OPE 
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FIG. 5: (color online) Spin observables for ^He two-body photodisintegration at photon laboratory energy — 12 MeV (left), 
— 20.5 MeV (middle) and E^ — 50 MeV (right). The upper rows show the analyzing powers for photon (^^(7)) and ^He 
(A^(^He)). The lower rows show spin correlation coefficients: Cxy and Cyx- The bands and lines have the same meaning as 
in Fig. H 

and short-range contributions not considered in the present work. We further emphasize that the results based on 
the AVI 8 potential and the corresponding MEG agree with the (present) chiral EFT calculation. 

For the three-body breakup of ^He, we only show the semiexclusive differential cross section d^cr/dQpdEp (where 
only one proton would be detected at 15 degrees with respect to the photon beam) at three photon laboratory energies 
Ej= 12, 20.5 and 50 MeV. The calculated cross section is shown as a function of the proton energies in Fig. [6l For 
the lower photon energy (left panel), the obtained bands appear to be relatively narrow, especially for higher proton 
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FIG. 6: Differential cross section for semiexclusive "^He three-body photodisintegration '^He(7,p)pn for proton emissions at 



— 15^ and photon laboratory energy — 12 MeV (left), 
and lines have the same meaning as in Fig. (4] 



^ 20.5 MeV (middle) and E^ = 50 MeV (right). The bands 



energies, where they both coincide with the AV18 results. For the lower proton energy the bands become broader. 
The TPE contributions bring the results close to the one obtained within the conventional framework. The situation 
is quite different for the highest photon energy (right panel). The shape of the calculated cross section is much more 
complicated in this case. Further, the band resulting from the TPE parts of the current operator appears to be very 
broad in the whole range of the proton energies. It remains to be seen whether the inclusion of the missing meson 
exchange current contributions will allow to reduce the theoretical uncertainty for the cross section. 



VI. CONCLUSIONS 



In this work we explored the effects of the TPE currents derived recently in the framework of ChEFT [16] in 
the deuteron and ^He photodisintegration reactions. We studied the role of various ingredients of the chiral 2N 
current operator in the unpolarized cross section and several polarization observables. As a main outcome of our 
study, we found that the new terms in the exchange current operator beyond the well-known one-pion exchange 
contribution play an important role for nearly all considered reactions. In particular, the differential cross section and 
the photon analyzing power in the deuteron photodisintegration process and the spin observables in ^He two-body 
photodisintegration are found to provide excellent testing ground for probing the fine details of the exchange current 
operator. 

We also found that the inclusion of the TPE contribution to the current operator alone typically results in very 
broad bands for the considered observables. This behavior is not unexpected. The OPE and TPE MEG are computed 
in the framework of chiral EFT within the low-momentum expansion and thus feature singular behavior at short 
distances (or large momenta). This leads to the observed large sensitivity of the calculated nuclear matrix elements 
to the short-distance behavior of the corresponding wave functions which is strongly scheme and cut-off dependent. 
In a complete calculation, the cut-off dependence of the low-energy observables is expected to be strongly reduced 
by the "running" of the corresponding short-range current operators, see e.g. Ref. [27] for the explicit examples of 
such a behavior in the case of the Ml properties of light nuclei within the hybrid approach and Ref. [28| for an 
extensive discussion on the (meaning of) renormalization in the context of nuclear EFT with a finite cut-off. Thus, 
the strong cut-off dependence in the obtained incomplete results which do not include the short-range contributions to 
the current operator should not be surprising. We expect that a complete NLO calculation including the short-range 
contact and the subleading OPE contributions to the 2N current operator will yield much narrower bands allowing 
for a quantitative description of electromagnetic reactions in a wider kinematical range. Work along these lines is in 
progress. 
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